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Abstract 

We consider the Laplace-Beltrarni operator in tubular neighbourhoods of 
curves on two-dimensional Riemannian manifolds, subject to non-Hcrmit- 
ian parity and time preserving boundary conditions. We are interested in 
the interplay between the geometry and spectrum. After introducing a 
suitable Hilbert space framework in the general situation, which enables 
us to realize the Laplacc-Bcltrami operator as an m-sectorial operator, 
we focus on solvable models defined on manifolds of constant curvature. 
In some situations, notably for non-Hermitian Robin-type boundary con- 
ditions, we are able to prove either the reality of the spectrum or the 
existence of complex conjugate pairs of eigenvalues, and establish similar- 
ity of the non-Hermitian m-sectorial operators to normal or self-adjoint 
operators. The study is illustrated by numerical computations. 
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1 Introduction 



Many systems in Nature can be under first approximation described by linear 
second order differential equations, such as the wave, heat or Schrodinger equa- 
tion. The common denominator of them is the Helmholtz equation describing 
the stationary regime and leading to the spectral study of the Laplace operator. 
Already from the mathematical point of view, it is important to understand 
the influence of the geometry to the spectrum of the Laplacian, subject to vari- 
ous types of boundary conditions, and vice versa, to characterize geometric and 
boundary interface properties from given spectral data. 

In this paper, we are interested in the interplay between the curvature of 
the ambient space and the spectrum of the Laplacian subjected to a special 
class of non-Hermitian boundary conditions. We choose probably the simplest 
non-trivial model, i.e., the spectral problem corresponding to the equation 

-A^ = X^ in n, (1.1) 

where A is a spectral parameter, fl is an a-tubular neighbourhood of a closed 
curve r (usually a geodesic) in a two-dimensional Riemannian manifold A (not 
necessarily embedded in M'^), i.e., 

n := {x e A \ dist{x,T) < a} , (1.2) 

and —A is the associated Laplace-Beltrami operator. The boundary conditions 
we consider are general 'parity and time preserving' boundary conditions intro- 
duced in Section |2.3.2| below; a special example is given by the non-Hermitian 
Robin-type boundary conditions 

^+iaip^O on dn , (1.3) 
on 

where n is the curve normal translated by geodesies orthogonal to F and a is a 
real- valued function. 

The Schrodinger equation in tubular neighbourhoods of submanifolds of 
curved Riemannian manifolds has been extensively studied in the context of 
quantum waveguides and molecular dynamics (c/. (5S] for a recent mathemat- 
ical paper with many references). Here the confinement to a vicinity of the 
submanifold is usually modelled by constraining potentials [40l [55] or Dirichlet 
boundary conditions [TH [331 [31] • 

Note that, on the contrary, the non-Hermitian nature of boundary condi- 



tions (1.3) enables one to model a leak/supply of energy from/into the sub- 
system f2, since the probability current does not vanishes on the boundary 90 
unless a = 0. In fact, non-Hermitian boundary conditions of the type ( |1.3[ ) has 
been considered in [30l [29l [28] to model open (dissipative) quantum systems. 
One also arrives at ( |1.3[ ) when transforming a scattering problem to a (non- 
linear) spectral one [T71 Ex. 9.2.4]. Finally, let us observe that Robin boundary 
conditions are known under the term impedance boundary conditions in classi- 
cal electromagnetism, where they are conventionally used to approximate very 
thin layers [HI [HI [6]. 
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Our primary motivation to consider the spectral problem ( |1.3[ ) comes 

from the so-called 'PT-symmetric quantum mechanics' originated by the pa- 
per [8] , where the authors discussed a class of Schrodinger operators H in (K) 
whose spectrum is real in spite of the fact that their potentials are complex. 
They argued that the rather surprising reality of the spectrum follows from the 
PT-symmetry property: 

[H,VT]^0. (1.4) 

Here the 'parity' P and 'time reversal' T operators are defined by {V'4)){x) := 
tp[—x) and Tip '■= ip- It is important to emphasize that T is an antilinear 



operator and that (1.4) is neither suflacient nor necessary condition to ensure 
the reality of the spectrum of H. 

Nevertheless, later on it was observed in [SI UTJ |33] that if the spectrum 
of a T^T-symmetric operator iJ in a Hilbert space T-L is indeed real (and some 



further hypotheses are satisfied) the condition (1.4 1 actually implies that H is 
'quasi- Hermitian' |50j . i.e., there exists a bounded invertible positive operator O 
with bounded inverse, called 'metric', such that 

H*^e-^HQ. (1.5) 

In other words, H is similar to a self-adjoint operator for which a conventional 
quantum-mechanical interpretation makes sense. We refer to recent reviews [71 
|44] and proceedings [27l \5\ [20] for further information and references about the 
concept of PT-symmetry. 

In addition to the potential quantum-mechanical interpretation, we would 
like to mention the relevance of T^T-symmetric operators in view of their recent 
study in the context of superconductivity [HlllHI, electromagnetism and 
fluid dynamics [Ml [111 [13 [l2] . 

A suitable mathematical framework to analyse T^T-symmetric Hamiltonians 
is either the theory of self-adjoint operators in Krein spaces [3SJl26j or the J-self- 
adjointness [TUj. The latter means that there exists an antilinear involution J 
such that 

H*=JHJ. (1.6) 



The concept (1.6) is not restricted to functional Hilbert spaces and it turns out 



that the majority of T^T-symmetric Hamiltonians existing in the literature are 



indeed J-self-adjoint. In general, however, the properties (1.4), (1.5) and (1.61 
are all unrelated [HIl [S3| . 



Summing up, given a non-Hermitian operator H satisfying (1.4), two fun- 
damental questions arises. First, 

1. is the spectrum of H real? 

Second, if the answer to the previous question is positive. 



2. does there exist a metric Q satisfying (1.5) 



It turns out that the questions constitute a difficult problem in the theory of 
non-self-adjoint operators. 
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For this reason, one of the present authors and his coauthors proposed in 
(see also |35j l an elementary one-dimensional T^T-symmetric Hamiltonian, for 
which the spectrum and metric are explicitly computable. The simplicity of 
the Hamiltonian consists in the fact that it acts as the Hamiltonian of a free 
particle in a box and the non-Hermitian interaction is introduced via the Robin- 
type boundary conditions ( |1.3 1 only. The model was later generalized to a 
two-dimensional waveguide in [lOj , where the variable coupling in the boundary 
conditions is responsible for existence of real (or complex conjugate pairs of) 
eigenvalues outside the essential spectrum (see also [37]). 

In this paper we continue the generalization of the models of [351 dD] to 
curved Riemannian manifolds. This leads to a new, large class of T'T'-symmetric 
Hamiltonians. Our main goal is to study the effect of curvature on the spectrum, 
namely the existence/absence of non-real eigenvalues and the metric. 

The organization of this paper is as follows. 

In the following Section [2| we introduce our model in a full generality, in the 
sense that the ambient geometry and boundary interaction of the spectral prob- 
lem (1.1) are described by quite arbitrary (non-constant and non-symmetric) 
functions. Our main strategy to deal with the curved geometry is based on the 
usage of Fermi coordinates. 

In Section [3j we use the framework of sesquilinear forms to define the 
Laplace-Beltrami operator appearing in (1.1) as a (closed) m-sectorial operator 
in the Hilbert space LP'{Vl). We also explicitly determine the operator domain 
if the assumptions about the geometry and boundary-coupling functions are 
naturally strengthen. Moreover, we find conditions about the geometry under 
which the operator becomes 'PT-symmetric (and T-self-adjoint). 

In order to study the effects of curvature on the spectrum, in Section [4] we 
focus on solvable models. Assuming that the curvature and boundary-coupling 
functions are constant, the eigenvalue problem can be reduced to the inves- 
tigation of (infinitely many) one-dimensional differential operators with PT- 
symmetric boundary conditions. Here the previous results [361 135] and the 
general theory of boundary conditions for differential operators [351 SS] are ap- 
propriate and helpful. In particular, since the P7~-symmetric boundary condi- 
tions are (except one case excluded here by assumption) strongly regular ones, 
it is possible to show that the studied one-dimcnsional operators are 'generi- 
cally' similar to self-adjoint or normal operators. However, it remains to decide 
whether this is true for their infinite sum, i.e., for the original two-dimensional 
Laplace-Beltrami operator. To answer this in affirmatively, it turns out that the 
J-self-adjoint formulation of PT-symmetry (c/. the text around (1.6)) is funda- 
mental, with J = T playing the role of antilinear involution. The properties of 
the solvable models are illustrated by a numerical analysis of their spectra. 

The paper is concluded by Section [5] where possible directions of the future 
research are mentioned. 
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2 Definition of the model 



We use the quantum-mechanical framework to describe our model. 



2.1 The configuration space 

We assume that the ambient space of a quantum particle is a connected complete 
two-dimensional Riemannian manifold A of class (not necessarily embedded 
in the Euclidean space M.^). Furthermore, we suppose that the Gauss curva- 
ture K of A is continuous, which holds under the additional assumption that A 
is either of class or it is embedded in M.^. 

On the manifold, we consider a C^-smooth unit-speed embedded curve F : 
Z, — > A^ with I > 0. Since F is parameterized by arc length, the derivative 
T F is the unit tangent vector of F. Let N be the unit normal vector of F 
which is uniquely determined as the C^-smooth mapping from [—1,1] to the 
tangent bundle of A by requiring that N{s) is orthogonal to T(s) and that 
{T{s),N{s)} is positively oriented for all s G [-1,1] {cf. [51 Sec. 7.B]). We 
denote by k the corresponding curvature of F defined by the Frenet formula 
S/tT = kN, where V stands for the covariant derivative in A. We note that 
the sign of k is uniquely determined up to the re-parametrization s i— > — s of the 
curve F and that k coincides with the geodesic curvature of F if ^ is embedded 
in R^. 

The feature of our model is that the particle is assumed to be 'confined' to 
an a-tubular neighbourhood of F, with a > 0. fl can be visualized as the set 
of points q in A for which there exists a geodesic of length less than a from q 
meeting F orthogonally. More precisely, we introduce a mapping C from the 
rectangle 

no -.^ (-1,1) X {-a,a) = Ji X J2 (2.1) 
(considered as a subset of the tangent bundle of A) to the manifold A by setting 

C{xi,X2) := expr(^j)(A^(a;i)a;2) , (2.2) 

where exp^ is the exponential map of A at q £ A, and define 

n:=C{nQ). (2.3) 

Note that Xi >-> C{xi,X2) traces the curves parallel to F at a fixed distance [x2[, 
while the curve X2 1— >■ C{xi,X2) is a geodesic orthogonal to F for any fixed xi. 
See Figure [TJ 



2.2 The Fermi coordinates 

Throughout the paper we make the hypothesis that 



£ : ilo — > is a diffcomorphism. 



(2.4) 



Since F is compact, (2.4) can always be achieved for sufficiently small a [cf. [23l 
Sec. 3.1]). Consequently, C induces a Riemannian metric G on fio, and we 
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Figure 1: Strip and boundary conditions 



can identiiy the tubular neighbourhood ft C A with the Riemannian manifold 
(Sloj G). In other words, O can be conveniently parameterized via the (Fermi or 



geodesic parallel) "coordinates" {xi,X2) determined by (2.2). We refer to [23l 
Sec. 2] and [24] for the notion and properties of Fermi coordinates. In particular, 
it follows by the generalized Gauss lemma that the metric acquires the diagonal 
form: 

where / is continuous and has continuous partial derivatives 82/, 9|/ satisfying 
the Jacobi equation 



/(•,0) = 1 
d2f{;0) = -K. 



dlf + Kf^O with { ^ ,)\[ ' (2.6) 



Here K is considered as a function of the Fermi coordinates {xi,X2). 
2.3 The Hamiltonian 

We identify the Hamiltonian H of the quantum particle in fl with the Laplace- 
Beltrami operator — in the Riemannian manifold (fig, G), subject to a special 
class of non-self-adjoint boundary conditions. 

2.3.1 The action of the Hamiltonian 

Denoting by G"-' the coefficients of the inverse metric and \G\ := det(G), 
we have 

-Ag = -\G\-^/^d,\G\^^^G''dj = -r^dif-^di - f-^d2fd2 . (2.7) 

Here the first equality (in which the Einstein summation convention is assumed) 
is a general formula for the Laplace-Beltrami operator —Aq expressed in local 
coordinates in a Riemannian manifold equipped with a metric G. The second 



equality uses the special form (2.5), for which |G| = and G~^ — diag(/~^, 1). 

Henceforth we assume that the Jacobian of ( |2.4[ ) is uniformly positive and 
bounded, i.e., 

/,/-! GL-(f^o), (2.8) 
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so that — Ag is a uniformly elliptic operator. Again, (2.8) can be achieved for 



sufficiently small a, cf. (2.6). 



Remark 2.1. The assumption (2.4) is not really essential. Indeed, abandoning 



the geometrical interpretation of 17 as a tubular neighbourhood embedded in A, 



{Qq,G) with (2.5) can be considered as an abstract Riemannian manifold for 
which (2.8) is the only important hypothesis. The results of this paper extend 



automatically to this more general situation. 
2.3.2 The boundary conditions 

We denote diflo — r,^ur^ the boundary in Xi direction, i e {1, 2}, see Figure[l] 



rf := {±1} X J2, rf Ji X {±a}. 



(2.9) 



Boundary conditions imposed respectively on diflg and ^2^0 ^^'c of different 
nature. Having in mind the situation when F is a closed curve, standard periodic 
boundary conditions are imposed on Oi^Iq, i.e., 



i>{-l, X2) = il^il, X2), diip{-l, xi) = ditpil, X2), 



(2.10) 



for a.e. X2 G t/2, where denotes any function from the domain of H. We 
assume also the symmetry condition on the geometry 



y{xi,X2) e : f{-l,X2) = f{l,X2), 



(2.11) 



(2.12) 



in order to have indeed periodic system in xi direction. 

On the other hand, non-self-adjoint PT-symmetric boundary conditions are 
imposed on 92 ^^o- A general form of PT-symmetric boundary conditions was 
presented in [5]; further study and more general approach to extensions can be 
found in [U |3] . Denoting 

^d2^) ' 

there are two types of the conditions, separated and connected. 
I. separated: 

(±P{xi) + ia{xi) 0^ 

for a.e. xi £ Ji, with a,/3 being real- valued functions. 
II. connected: 



*(xi,±a) = 



(2.131) 



*(xi,a) = B[xi)-^{xi,-a), 
for a.e. xi € Ji, where the matrix B has the form 



(2.13 II) 



B{xi) 



^1 + 6(a;i)c(xi)e''^(^i) b{xi) 

c{xi) + K2;i)c(xi)e-'"*(^i) 
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with b,c,<j> being real- valued functions satisfying b > 0, c > —1/b, (j) G 
[-7r,7r). 

We specify assumptions on smoothness, boundedness and periodicity of the 
functions entering the boundary conditions later. The index t e {1, 11} wiU be 
used throughout the paper to distinguish between the two types of boundary 
conditions. 

The boundary conditions (2.13.) are T'T-symmetric in following sense: if a 
function t/j satisfies (2.13 ,), then the function VTip satisfies ( 2.13[ .) as well. Here 
and in the sequel the symmetry operators V and T are defined as follows: 



(2.14) 



It is important to stress that the T'T-symmetric boundary conditions ( 2.13[ .) 
do not automatically imply that the operator H is T^T-symmetric, unless addi- 
tional assumption on the geometry of i7o is imposed. The assumption, ensuring 
the T'T-symmetry of H {cf. Proposition 3.1 below), reads 



V(a;i, a;2) G : f{xi,X2) = f{xi,~X2). 



(2.15) 



In view of (2.6 1, a necessary condition to satisfy the second requirement in (2.151 
is that the curve F is a geodesic, i.e. k = 0. 



2.3.3 The functional spaces 

The space in which we give a precise meaning of H is the Hilbert space L^(fio: G), 
i.e., the class of all measurable functions ip, ip on J7o for which the norm || • \\q 
induced by the inner product 



is finite. Assuming (2.8 1, the norm 

in L^(rio)- Moreover, the 'energy space 



f{x)^p{x) \G{x)\^'^dx (2.16) 
• II G in i^(rio, G) is equivalent to the usual 



one 



M^i'2(17o,G) := {V e £2(r!o,G) | |VgVIg - 9,VG'^a,V e L\no,G)] (2.17) 

can be as a vector space identified with the usual Sobolev space 14^^'^ (rig)- 

However, this equivalence does not hold for M^^'^-spaces, unless one assumes 
extra regularity condition on /: 



yx2 £ J2 ■■ f{;X2), .rH;X2) e Tyi'°°(Ji). 



(2.18) 



Under this assumption, which is actually equivalent to the Lipschitz continuity 
of /, /^^ in the first argument {cf. [111 Chapt. 5.8.2.b., Thm. 4]), one can indeed 
identify the ly^'^-Sobolev space on the Riemannian manifold {^Iq, G) (precisely 
defined, e.g., in [25, Sec. 2.2]) with the usual Sobolev space W'^''^{ilo). 
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2.3.4 The schism: two definitions of the Hamiltonian 



Althougti the above equivalence of the M^ ' -spaces under the condition (2.181 



is not explicitly used in this paper, it is in fact hidden in our proof that the 
particle Hamiltonian on L^(Oo,G') naturally identified with 

H,^b := -AgV", (2.19a) 
iP e Dom(iJ,,) := {V' e W^-^irio) \ ijj satisfies plof and (|2.13[,)}. (2.19b) 



is well defined (c/. Theorem 3.1 1. As mentioned in Section 2.3.2[ we use the 
notation i/^, with t £ {I, H}, to distinguish between separated (2.131) and 



connected (2.13 II) boundary conditions. 



To avoid the additional assumption (2.181, one can always interpret (2.7) in 



the weak sense of qua drat ic forms, which gives rise to an alternative Hamilto- 
nian (c/. Corollary 3.1). This is the content of the following section, where 
we also show that H^^ ^ provided that (2.111, (2.181, and some analogous 



hypotheses about the boundary-coupling functions hold. 



3 General properties 



The main goal of this section is to show that the Hamiltonian introduced 



in (2.19) is a well defined operator, in particular that it is closed. This will be 



done by proving that — H,^, where is the alternative operator defined 
through a closed quadratic form. Finally, we establish some general spectral 
properties of the Hamiltonians. 



3.1 The Hamiltonian defined via quadratic form 

Taking the sesquilinear form {lp^H^iI))q with ip,ip £ Doni(iJt) and integrat- 
ing by parts, one arrives to a sesquilinear form, which is well defined for a 
wider class of functions (p,ip, not necessarily possessing second (weak) deriva- 



tives. The function / is assumed to satisfy (2.81 and (2.11), however the extra 



regularity condition (2.18) is not required 



More precisely, exclusively under assumption (2 
the sesquilinear form 



for a moment, we define 



where, for any Lp,ip £ Dom (/i,,). 



(52<p,52'0) 



G ' 



((^, (/3 -t- ia)ijfj + ((^, (/3 - ia)ij)\ 
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Here Bij denotes the elements of the matrix B defined in ( 2.13 1, the operator V 



is introduced in (2.14) and 



iip,ip)c J ipixi,±a)ip{xi,±a) f{xi,±a)dxi. 
All the boundary terms should be understood in sense of traces [T]. 



Lemma 3.1. Let f satisfy (2.8). The forms h^^h^ are densely defined, is a 



symmetric, positive, closed form (associated to the self-adjoint Laplace- Beltrami 
operator in L^iyi^, G) with periodic boundary conditions on diVl^) and Neumann 
boundary conditions on d2^t))- 

Proof. The density of the domains is obvious, properties of are well known, 
see the detailed discussion on a similar problem in [15l Sect. 7.2]. □ 

Although the forms are not symmetric, we show that can be understood 
as small perturbations of . 



Lemma 3.2. Let b,l/b,c,a, j3 € L°°{Ji) and let f satisfy (2 
relatively bounded with respect to with 



\hl[n<eh^[^l,]+e-^CU\\ 



2 



Then h, are 



(3.1) 



for all -0 e WjJ^J(r2o) and any positive number e. The constant C depends on 
the function f , dimensions a, I, and boundary- coupling functions a, /3 or b, c, (j). 

Proof. The proof is based on the estimate 



(3.2) 



where e is an arbitrary positive constant and C is a positive constant depend- 
ing only on a and I. We give the proof for hf only because the other case is 
analogous. The assumptions on a, /3 and property (2.8 1 allow us to estimate 
the functions \a\, \f3\ and / by their i°°-norms. Consequent application of (3.2 1 
therefore yields 

In order to replace the term jjV'i/'lP by h^[4'], the regularity assumption on 

and 

□ 



geometry (2.8) is used. Once we consider the equivalence of the norms 



|g and the arbitrariness of e, we obtain the estimate (3.1) 



Corollary 3.1. Let 6, 1/6, c, a, /? e L°°{Ji) and let f satisfy (2.8). Then there 
exist the unique m-sectorial operators in L'^{Q.q,G) such that 

K{ip,i;)^:{ip,H,i;)G (3.3) 
for all Tp G Dom (iJj and ip G Dom {h^), where 

Dom(4) {0 e W^HQo) \ 3F e L\no,G), e W^iXQa), 

Ki^,i;)^{cp,F)G}. (3.4) 
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Proof. With regard to Lemmata 3.1 3.2 and the perturbation result [5T1 Thm. VI. 3. 4] . 
the statement follows by the first representation theorem fSP, Thm. VI. 2.1]. □ 



3.2 The equivalence of the two definitions 

Under stronger assumptions on smoothness of functions appearing in boundary 
conditions (2.13 .) and on the function / entering the metric tensor G, we show 
that operators i?^ associated to the forms are equal to the Hamiltonians H,^ 
defined in ( 2. 19 ) . To prove this, we need the following lemma. Let us introduce a 
space of Lipschitz continuous functions over [— Z, I] satisfying periodic boundary 
conditions: 



l,oo 



per 



(Ji) :={^eiyi'°°(Ji)| ^(-0 = ^(0}. 



Lem ma 3 .3. Let a, I3,h,\/h,c,4> £ M^pc?°(^i) and let f satisfy ( |2.11[ ), 
and (2.18). Then for every F G L'^{flQ,G), a solution Tp to the problem 



belongs to Dom(iJt) introduced in ( 2.19b[ ). 



(3.5) 



Proof. We prove the separated boundary conditions case only, the connected 
case is analogous. For each ip e Wp^{no) We introduce a difference quotient 



(3.6) 



where ips{xi,X2) := i^{xi + i5, ^2) and (5 is a small real number. The shifted 
value -015(3^1; 3^2) is well defined for every xi G Ji and 5 G K by extending ip 
periodically to R. We use periodic extensions of other functions in Xi direction 
throughout the whole proof without further specific comments. The estimate 



(3.7) 



is valid for 6 small enough [19, Sec. 5.8.2., Thm. 3]. 

We express the difference of identities (3.5 1 for ijj and ipg, whence we get for 
every </7G VFpi'?(17o) 

{Sf-^)d,^j) + {d,^,fi'di{6tlj)) + {82^, {df)d2ij) 

+ {d2V, f&d2m) + + '^a))i^sf^' + /(/3 + i«)(5V))'^^ 

+ 5{f{P - iamf^ + (vp, /(/3 - ia){Si;)f^ 

^{ip,i6f)Fs) + {ipjm), (3.8) 
where (•, •) is the inner product in L'^lflo) and 



ip[xi,±a) ip{xi, ±a) dxi. 



(3.9) 
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We insert ip = 54> into equation (3.8 1 and apply the 'integration-by-parts 



formula [13 Sec. 5.8.2] for difference quotients, i.e., {ip,SF) = —(^{—6)(f,F), 
in order to avoid the difference quotient of the arbitrary {e.g. possibly non- 
continuous) function F E L^{nQ,G). Using the embedding of W'^''^{flo) in 
L^(9r2o)j the regularity assumptions on a, P and /, the Schwarz and Cauchy 
inequalities, and the estimate (|3.7|), we obtain 



(3.10) 



where C is a constant independent of 5. By standard arguments [191 D.4], this 
estimate yields that g M^"'^'^(ilo). 

At the same time, standard elliptic regularity theory [2^ Thm. 8.8] implies 



that the solution t/i to (3.5 1 belongs to W^^^ {^a). Thus ■0 satisfies the equation 

-AgV' = -F' (3.11) 



a.e. in JIq. If we express 91-0 from (3.11), we obtain that (9|?A G L?{Q,q). 

It remains to check boundary conditions of Dom(i7i). Once the W^'^- 
regularity of the solution -0 is established, this can be done by using integration 



by parts in the identity (3.5) and considering the arbitrariness of tp, see [101 
Lemma 3.2] for the more detailed discussion in an analogous situation. □ 

Let us write Hi{a, l3) and Hii(b, c, </>) if we want to stress the dependence of 
the Hamiltonians on functions a, /3 and 6, c, (j) entering the boundary conditions. 



The orem 3.1. Let a, ^, 6, 1/6, c, G Wli^{Ji) and let f satisfy ( |2.11 ) 



and (2.18). Then 

1. He = H„ 

2. are m-sectorial operators, 

3. the adjoint operators H* can be found as 



4-. the resolvents of are compact. 

Proof. Ad 1. It is easy to verify, by integration by parts, that if G Dom(iJt) 
then G E)om(iJt); in fact, the function F from (3.4) satisfies F — —Acip in 
the distributional sense. Thus i?,, C i/f The more non-trivial inclusion C 
follows from Lemma [3. 3[ Once the equality of the operators is established, the 



other properties readily follow from the corr esponding properties for H^. 



3.1 



Ad 2. is m-sectorial by Corollary 
Ad 3. By Thm. VI. 2. 5], the adjoint operator i/* is associated to the ad- 
joint form hl{ip, ■0) :— h^{ip, ip), which establishes the required identities for H^. 
Ad 4- The compactness of the resolvents for is provided by the pertur- 



bation result |3]j Thm. VI. 3. 4] and Lemmata 3.1 3.2 



□ 
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3.3 Spectral consequences 



Since the Hamiltonians are m-sectorial by Theorem 3.1 the spectrum (as 
a subset of the numerical range) is contained in a sector of the complex plane, 
i.e., there exists a vertex 7 G K and a semi-angle 9 S [0,7r/2) such that 

a(i/.)c {CeC I |arg(C-7)l <^}. 

Furthermore, since the resolvents of are compact, the spectra of iJ,, are purely 
discrete, as it is reasonable to expect for the Laplacian defined on a bounded 
manifold. 

Under the additional assumptions on the geometry of the model (2.15), one 
can show that arc PT-symmctric. 



Proposition 3.1. Let a, /3, &, 1/6, c, e W^^i^Ji) and let f satisfy (2.8), 
(2.15), and (2.18). Then Hamiltonians H,^ are 

1. VT -symmetric, i.e., {VT)H, C H,{VT), 

2. V -pseudo-Hermitian, i.e., = VH*!^, 

3. T -self-adjoint, i.e., = TH*T, 



where the operators V and T are defined in (2.14). 

Proof. Note that the T'T-symmetry relation means that whenever ip € Dom {H^), 
VTip also belongs to Dom{H,) and VTH^ij^ ~ H^VTip. This can be verified 
directly using the definition of iJ^ via (2.19). The proofs of the remaining 
statements are based on the explicit knowledge of the adjoint operators. Theo- 
rem 03. □ 



Corollary 3.2. Under the hypotheses of Proposition \37l\ the spectra of are 
invariant under complex conjugation, i.e., 



VA e 



A e (j{H,) 



Proof. Recall that the spectrum of is purely discrete due to Theorem 3.1 4. 
With regard to T'T-symmetry, it is easy to check that if ip is the eigenfunction 
corresponding to the eigenvalue A, then VTip is the eigenfunction corresponding 
to the eigenvalue A. □ 



4 Solvable models: constantly curved manifolds 

In order to examine basic effects of curvature on the spectrum of the Hamilto- 
nians we investigate solvable models now. We restrict ourselves to the spectral 
problem in constantly curved manifolds and subjected to constant interactions 
on the boundary, i.e., the functions K,a, /3,b,c,4> are assumed to be constant. 



Moreover, we assume that F is a geodesic, i.e. k — 0, to have (2.15) 
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4.1 Preliminaries 

To emphasize the dependence of the Hamiltonians on the curvature K, we use 
the notation H^(^x^ in this section. One can easily derive the scahng properties 
of eigenvalues for constant K ^ 0: 

K{K,a,l) = |i^|A,(±l,v^a,v^O- 

Hence, the decisive factor for qualitative properties of the spectrum is the sign 
of K, while the specific value of curvature is not essential. Hereafter we restrict 
ourselves to 

if e {-1,0,1}. (4.1) 

Possible realizations of the ambient manifolds A corresponding to these three 
cases are pseudosphere, cylinder, and sphere, respectively, see Figure |2] 




(a) K = —I, pseudosphere (b) K = 0, cylinder (c) K = 1, sphere 

Figure 2: Realizations of the constantly curved manifolds. 



Remark 4.1. The pseudosphere should be considered as a useful realization 
of A with K — —1 only locally, since no complete surface of constant negative 
curvature can be globally embedded in (this is reflected by the singular 
equator in Figure [2] (a)). However, since O is a precompact subset of A, the 
incompleteness of the pseudosphere is not a real obstacle here. 

Moreover, hereafter we put Z = tt, so that the length of the strip is 2tt. 
This provides an instructive visualization of f2 as a tubular neighbourhood of a 
geodesic circle on the cylinder and the sphere, see Figure [2j 

For K = and constant curvatures (4.1), the Jacobi equation (2.6) admits 
the explicit solutions 



f(K){xi,X2) 



cosha;2 if K = —1 . 
1 if if = , 

cos X2 if K = 1 . 



(4.2) 
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It follows that the assumption (2.8) is satisfied for any positive a if = —1, 0, 
while one has to restrict to a < tt/2 ii K = 1. The latter is also sufficient to 



satisfy (2.4) for the sphere. There is no restriction on a to have (2.4) if F is 



the geodesic circle on the cylinder. In any case (including the pseudosphere) 



(2.4) can be always satisfied for sufficiently small a. The other hypotheses, 



i.e. (2.11), (2.15), and (2.18), clearly hold regardless of the curvature sign. 



Remark 4.2. In view of Remark 2.1 a < 7r/2 for K = I is the only essential 
restriction in the constant-curvature case (4.2). 



Explicit structures of the Hamiltonians H,(^x) introduced in (2.19) readily 
follow from (2.7) by using (4.2): 



\ — ai^-9^-tanhx292 if K = -1 , 

cosh X2 



H,(K) = < -dl - d, 



if K = 0, 



(4.3) 



^ — df -di+ tan x^d^i if K = l, 

COS^ X2 



on V)ou\{H,(K)) 



4.2 Partial wave decomposition 

Since both the coefficients of ) and the boundary conditions are independent 
of the first variable xi, we can decompose the Hamiltonians into a direct sum 
of transverse one-dimensional operators. The decomposition is based on the 
following lemma. 



Lemma 4.1. 

V*e L2(f]o,G), 

where 



'i!(xi,X2) ^^^„^{x2)(j>v^{xl) in L'^{nQ,G), (4.4) 



mG2 



0m(a;i) 



V27r 



1pmix2) := (0m, *(-,.T2)) 



(4.5) 



Proof. We may restrict the proof to L^{fto) only because the norms || • || and 
II • IIg are equivalent due to (|2.8|. Let us also stress that G is independent of xi 
and {0m}^ 



forms an orthonormal basis of L (Ji). Hence 



1pm{x2)<l>r, 



\\^i;X2)\\LHJ,)eL^iJ2). 



(4.6) 



The decomposition in L^(ilo) can be then justified by using the dominated 
convergence theorem. □ 



16 



Writing 'i'{xi,X2) = Smsz </'™(2^i)V'm(2:2) in the expression H,(x)'^ and 
formahy interchanging the summation and differentiation, we (formaUy) arrive 
at the decomposition: 

H.iK) = ^r(V)5" (4.7) 



with 



H. 



,{K) - 



—82— tanhx292 H 2 — 

cosh X2 

— 9| + m? 

9 

— 9| + tanx292 



if K = -\, 
if K^O, 
if K^l, 



cos^ 0:2 

where B"^ are bounded rank-one operators defined by 

X2) := ((/),„, V'l', a;2))^2(j^) '/'m(a;i) • 
The operators -ff™^) act in J2, di'(i^)) spaces with the measure 

{coshx2da;2 if if = — 1 , 
da;2 if if = , 

cos a;2 dx2 if if = 1 . 

The domains of HJ^^^^ are given by 

Dom (ii7(V)) := l^/- e iy2'^(J2) I satisfies (|2.13[^)}, 



(4.8) 



(4.9) 



(4.10) 



with obvious modification of the T'T-symmctric boundary conditions ( 2.13[ .) 
to the one-dimensional situation. 

To justify the decomposition (4.7) in a resolvent sense, we need the following 
technical lemma specifying the numerical range of i?™;^) • 

Lemma 4.2. Let denote the numerical range oj H^^y Then for every 

m^O there exist real constants Co,Ci independent of m such that 



C |z e C I Rez > co + rn^, |Imz| < ci^jRez^ |co| - m^j . (4.11) 



Proof. We give the proof for ^^i^+i) "-"^1^' other cases are analogous. We 
abbreviate (•,•)+ := (■, ■)L2(j2,di.(+i)) and define 



V"'{X2) 



COS-^ X2 



h[i,] := (V,iii7+i)^)^ 



for every € Dom (iij"Yjj). Integration by parts yields the following expressions 
for real and imaginary parts of h[ip]: 

Reh[^] = ||V''||'+ + (V,«"V)++/?cosa(|V'(a)|2 + |V'(-a)P), 
lmh[ip] = acosa[\'ilj{a)f — \'ip{—a)f) , 



17 



for every ip S Dom {H^^^^). The estimates of Re h[ip] and Im h[ip] can be easily 
obtained taking into account the equivalence of the norm || • ||l2(J2) with || • 11 + 
and using the one-dimensional version of the estimate (3.2). □ 

Now we are in a position to establish the main result of this subsection. 
Proposition 4.1. Z) := ^?(^t(V)) non-empty and D C f?(i?t(A'))- For 



every z Cz D, 



{H,^K)-zr' = ^iH:i^i,^-z)-'B^' 



(4.12) 



where {H™j^^— z) ^ abbreviates l(E){H™j^^—z) ^ acting on L^{Ji)®L?{J2,Avi^j^)) 
and _B™ are defined in (4.8). 

Proof. We give a proof for H^j^^^ only, the remaining cases are analogous. Take 
z G D, for every G L^{^, G) and m € Z, we define 



Um{x2) 



(4.13) 



where ipm was introduced in (4.5). It is clear that Um G L^( J2, dz/(+i)). With 
regard to Lemma 4.2 take tuq G Z such that for every m > toq: z ^ "^"if)- 
Using [3TJ Thm. V .2.3] together with Lemma 4.2 we get for m > mo 

ll'0m||L2(J2) 



1 



(4.14) 



where Ci is a constant independent of m, nonetheless depending on z, |a|, |/3|, 
and a. Let us remark that since z G 13, ||C/m|lL2(j2) are bounded for finitely 
many m smaller than uiq. From the identity 

ll^^mlll + (-ia + /3)cosa|t/™(a)p + {-ia - (3) cosa|?7™(-a)p 
+(i;"[/„,C/™)+-z||(7™||^ = (V,n,C/,„)+, 



with and (•, •)+ defined in Lemma 4.2 we obtain the estimate for the norm 
of for m > mo, 



+ 1 



Again, for finitely many m < mo, ||^^mllL2(j2) are clearly bounded. With regard 
to (4.6), (4.14), and (4.15), every function Rm{xi,X2) ■= (t)m{xi)Um{x2) belongs 
to W^?(f^o). 

Our goal is to show that R := J2mez-^"i is in Wp^^{ilo) as well. 

The finite number of bounded terms with m < tuq is included in the following 
estimates and equalities without any other specific comments. The identity (4.6 1 
and inequality (4.14) together with Fubini's theorem imply 
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A similar estimate can be obtained for d2Rm provided that we use the inequal- 
ity (4.15). For diR,n, we have 

2 



N 



J2 ^1^" 



m=-N 



N 



N 



"l|lL2(.72)> 



m=-N 



m=-N 



where we used the inequality (4.14). The fraction in the sum on the right hand 
side is bounded, therefore, using the Parseval equality, the limit X^msz^i^™ 
remains in L^(r2o). We conclude that R belongs to VF^'^(rio) and 

\\R\\w^-''{no) ^ C'3||*||l2(o„). 

It remains to verify that R belongs to Wp^{VIq). We introduce the partial 
sum Rn ■■= YZ,=-nR^- The fact that Rjq e Wl^^in^) for every iV e N and 
the (trace) embedding of M^^'^(f2o) in L^(9rio) yields 

I (^, •) - R{1, •)) + 1 = I {v. Ri-i. •) - Rn{-1, •) + RnH, •) - •)) + 1 

< 2C4 \\R - i?Jv||wi-2(Oo) 

for every (p € L'^{J2,'^^(+i))i C4 is a constant depending only on V,q. Notice 
that the left hand side does not depend on N. Hence, the periodicity of R is 
justified by taking the limit N — +00 and considering the arbitrariness of (f. 
Now, knowing that R belongs to Wp^^{^lo), one can easily check that 



This implies that R E Dom see Lemma 

Proposition 4.1 



hi{ip, R) - z{ip, R)L^no,G) = iv, *)L2(no,G) • 

and (iJi(+i)-z)i? = □ 



3.3 



Corollary 4.1. 



has the important consequence for the spectrum of H^(^x)- 



Proof. The inclusion a(H^(^j^^^ C a(^Hl^j^-^^ follows from Proposition 
the other one is trivial. 



4.1 



Remark 4.3. Notice that the statement of Corollary |4. 1| relating the spectra of 
a direct sum of operators with their individual spectra does not hold in general 
(c/. [ini Thm. 8.1.12]). In our case, however, we have been able to prove the 
result due to the compactness of resolvents and additional information about 
the behaviour of the numerical ranges of HJ^^^-^ {cf. Lemma 4.2). 
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4.3 Similarity to self-adjoint or normal operators 

We proceed with an analysis of H^^^y For sake of simplicity, we drop the 
subscript 2 of the X2 variable in the sequel. We remark that PT-symmetry 
and T'-pseudo-Hermiticity of H^j^^ is preserved with V and T being naturally 
restricted to L-^( J2, di/(/f)). 

The operators -H^™^) neither self-adjoint nor normal, nevertheless we can 
show the following general result: 

Theorem 4.1. For every m E Z and K G { — 1,0, 1}; 

1. The families of operators H^j^-^{a, f3), i?^^^^ (6, c, 0) are holomorphic with 
respect to parameters a, (3, ana b, c, (j) entering the boundary conditions. 

2. The spectrum of H^-j^^ is discrete consisting of simple eigenvalues (i.e., 
the algebraic multiplicity being one), except of finitely many eigenvalues of 
algebraic multiplicity two and geometric multiplicity one that can appear 
for particular values of a, /3 and b, c, 4>. 

3. If all the eigenvalues are simple, then 

a) the eigenvectors of ^[(V) /c"^ a Riesz basis, 

b) H'^j^^ is similar to a normal operator, i.e., there exists a bounded 
operator g with bounded inverse such that gHJ^j^-^g^^ is normal, 

c) if moreover all eigenvalues are real, then H^j^^ is similar to a self- 
adjoint operator, i.e., gHJ^j^.^g~^ is self-adjoint. 

4. Let us denote by {V-'i.m}jgfj the eigenf unctions of H^j^y The set of eigen- 



functions :~ {0m''/'i.m}„g2 ieN' '^here <j)„i were introduced in (4.5) 
forms a Riesz basis of L^{Q,o,G). 

Proof Ad 1. In view of [SH Sect. VII, Ex. 1.15], the Hamiltonians iJ™^)(a,/3), 
considered as a family of operators depending on parameters a, /? entering 
boundary conditions, are holomorphic. The same is true for H^^^j^-^{b,c,(j)). 

Ad 2. The separated boundary conditions belong to the class of strongly 
regular boundary conditions 45, 4(1. The connected T'T-symmetric boundary 
conditions are strongly regular as well because 9i = —b, 6'_i = b (in Naimark's 



notation) and b is non-zero by the assumption in (2.13). Moreover, all the 
eigenvalues are simple |39j up to finitely many degeneracies that can appear: 
eigenvalues with algebraic multiplicity two and geometric multiplicity one. 

Ad 3. With regard to the strong regularity of boundary conditions, the 
eigenfunctions of the Hamiltonian H^^x) fo^'™ a Riesz basis [35], except the 
situations when the degeneracies appear. The existence of Riesz basis implies 
the similarity to a normal operator and as a special case the similarity to a 
sclf-adjoint operator if the spectrum of HJ^k) ^'^sd. 

In more details, let {V'nj^gpj be the Riesz basis of eigenvectors of W^j^y 
i.e., H™j^^ipn — ^nipn- By definition, there exists a bounded operator p with 
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bounded inverse such that {pV^nl^gp^ is an orthonormal basis that we denote 
by {en}„eN- Then 

is a normal (self-adjoint if every A„ € M) operator. 

Ad 4. At first we show that ^ is complete, i.e., SS^ = {0}. Take w G 
i.e., for every to G Z, i e N, 

0m(a;i)V'i,m(a;2)w(a;i,a;2)dxidi/(if)(x2) 



l/'»,mWm(x2)dl/(;f)(x2), 

where a;m(a;2) := / (/)m(a;i)aj(a;i, X2)da;i. Since {V'i,m}jgpj forms a Riesz basis, 
.J J\ 

uj,n — a.e. in J2, di^(x)) for every 171^X1. Since {0m}„jgz is the orthonor- 
mal basis of Ji), w = in L^(rio, G). 

Now we define an involution {ViiIj){xi,X2) '■= X2). We show that 

V'i.TO can be normalized in such way that ^ is T'lT-orthonormal, i.e., 

Since ViT4>m — 4>m, T^iT-orthogonality follows immediately for m ^ n because 
4>m are orthonormal in (•, •)l2(Ji) and G is independent of xi. For m = n we 
have 

{4'm'4'i,m,'PlT(f>m4'j,m) j^2i^^g Q) = (''/'i,'" J '^V'j,m) ^2 ( ,di/(ff ) ) " (4.16) 



If i 7^ J, then the right hand side of (4.16 ) is zero because T'4'j,m is an eigenfunc- 
tion of (^H™^jY . Indeed, it is a general fact that eigenfunctions of H and H* 
corresponding to different eigenvalues are orthogonal. It remains to verify that 
if i = j, then the right hand side of (4.16) does not vanish, i.e.. 

However, this is precisely the condition on Xj.m being a simple eigenvalue of 
^^K) ■ '^^^ b® easily seen either directly or it follows from [3T1 Thm. 5] . □ 

Remark 4.4. Notice that an additional symmetry with respect to Vi was es- 
sential in the proof. The set of eigenfunctions is not T-orthonormal because 
the products {(f>m''Pi,m,T4'm'4'i,m) j^2^Qg q-^ vanish. This situation is typical for 
T-self-adjoint operators with eigenvalues that are not simple j2T| . 
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4.4 Separated boundary conditions 

At first, we investigate the Hamiltonians i/j"Q-,(a, /3). Then H^_^^Ja,0) are 
analysed. These resuhs together aUow us to describe the remaining p ^ case. 

4.4.1 Zero curvature 

As expected, the zero curvature case is the simplest and it will serve as a refer- 
ence model. In fact, the corresponding one-dimensional eigenvalue problem 

f --0" + = fc^V" in (-a, a), 

has been already studied previously in |36j . Here we overtake the main results. 

Proposition 4.2. The spectrum of H"^Q-^{a,0) is real for all m E Z. The 

eigenvalues Xj^m and eigenf unctions 'il'j_m can he written in the following form, 
TO G Z, 



A.™ - " T"; V.::' (4.18) 




{Co exp (-iax) if j = , 

C. cos(k.x) + 77^-^ -. . , sm(fcox) if j > I, 

where kj '■— If oP' ^ k'^ , i.e., there is no level-crossing for the same to, 
then the operator is similar to a self-adjoint operator or, equivalently, it is 
quasi- Hermitian. 

Remark 4.5. Closed formulae for the metric operator Q for H^g^{a, 0) are pre- 
sented in |361 135j . The similarity transformation g can be found as g = \/Q or 
as any other decomposition of the positive operator Q — g*g. 

The a-dependence of eigenvalues A for to = 0, 1, 2 is plotted in Figure 3] 
The case of /3 7^ is more complicated and as it was remarked in [36], the 

spectrum of H^|^^{a,|3) can be complex. More precise results follow from a 

further analysis, not presented in [55] . 

Proposition 4.3. 

1. If P > 0, then the spectrum of II^Q^{a, {3) is purely real for all m E Z and 
aeR. 

2. If P < 0, then the spectrum of H^Q^{a, (5) is either purely real or there 
is one pair of complex conjugated eigenvalues with real part located in 
the neighbourhood of cP' -V fP' ■ For fixed negative j3, the points an where 
a pair of complex conjugated eigenvalues appears (by increasing of a) are 
determined by Q^-^-\~j3'^ — k'^, where fc^ := ((2rt + l)7r/4a)^ for some rt e N. 
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Figure 3: a-dependence of eigenvalues, zero curvature, a — 7r/4. Red, green, 
and blue colour corresponds to m = 0, 1, 2 respectively. 



The eigenvalues X = of H^^^^^{a, j3) are determined (k — is admissible only 
if a 



= 0) by the equation 

{k^ ~ sin(2fca) - 2pk cos(2/ca) = 0. 
The corresponding eigenf unctions read 



ip{x) — C cos{kx) + 



k sin(fca) — (ia + f3) cos(fca) . 



k cos(fca) + (ia + f3) sin(fca) 



sin(A:a;) 



(4.19) 



(4.20) 



The eigenvalues of H^^^ (a, /3) are obtained by adding to the eigenvalues of 
7?i%)(a,/3). 

Proof. We proceed in a similar way as in the alternative proof [361 Sect. 6.1] 
of the reality of the spectrum of H^Qj{a,0). The original eigenvalue problem 



(4.171 with m = can be transformed, using 4'{x) := e'"^7/)(a;), into 



— 0" + 2ia0' + = A(/) in (—a, a), 
0'(±a) ± /?0(±a) = 0. 



(4.21) 



We multiply the equation (4.21) by 0" and integrate over {—a, a). Next we 
multiply the complex conjugated version of the equation (4.21 ) by 0" and again 



we integrate over (—a, a). By subtracting the results and integrating by parts 



with use of the boundary conditions in (4.21), we obtain the identity 

- a/3' (|0(a)p - |0(-a)p) = Im A + /?(|0(a)p + \^{~a)\') 

(4.22) 

If we perform the same procedure, however, with multiplication by (f>, after some 
integration by parts we receive the relation 



(|0(a)|2-|^(-a)P) =ImA||0||i.(,,^). 



(4.23) 



23 



Combining (4.22) with (4.231 leads to the identity 



= IniA 



(4.24) 



If /3 is positive, then the whole term in the brackets is strictly positive and thus 
imaginary part of A must be zero. This proves the first item of the proposition. 
If /3 is negative, then complex eigenvalues can appear. If we divide the 



equation (4.19) by fc^ and leave only sin(2fca) term on the left hand side, then it 
is clear that eigenvalues approach {rm/iaf' for k real and large enough. After 



simple algebraic manipulation (4.19) becomes 

213k 



tan(2fca) 



(4.25) 



fc2 - a2 _ p2 

and the eigenvalues correspond to the intersections of the graphs of functions on 



left and right hand side of (4.25). We denote l{k) the function on the left hand 
side, r(fc) the one on the right hand side, and fco :— \Ja^ -\- fi'^. The behaviour 
of r(fc) for /j e M is summarized in Table [T] 



k 


sign 


asymptotics 


(-0O, -fco) 


r{k) > 0, /(/c) > 0, > 


limfe^_fe„_ r{k) = +oo 
limfe^_oo r{k) 


(-A:o,0) 


r{k) < 0, r'(fc) > 0, r"{k) < 


\h-nk^-ko+ r{k) = -oo 


(O,fco) 


r{k) > 0, r'(fc) > 0, r"(fc) > 


limfe^fe„_ r{k) = +oo 


(fco,oo) 


r{k) < 0, r'(fc) > 0, r"(fc) < 


liuik^ko^ r{k) = -oo 
limfe^oo r{k) = 



Table 1: The behaviour of r{k). 



Graphs of functions l{k) and r{k) are plotted in Figure [i] It is clear from the 
holomorphic dependence of eigenvalues on a, ^ (a consequence of Theorem |4.1[ ) 
that eigenvalues are close to {mr/a)'^, corresponding to zeros of l{k) — tan(2fca), 
except those in the neighbourhood of + /J^ . Since a = (3 = case corresponds 
to Neumann boundary conditions, for small a and /?, all eigenvalues must be 
close to (n7r/2a)2. Hence, if we fix /? and increase a, then two intersections of 
graphs of l{k) and r(fc) are "lost" precisely at the point where + = k'^ 
for some n € N, i.e., the asymptote of r(/c) corresponds to the asymptote 
of the tangent l{k). This implies the creation of complex conjugate pair of 
eigenvalues. If we increase a more, two intersections appear again which means 
the annihilation of complex conjugate pair, i.e., the restoration of two real 
eigenvalues. The two intersections are lost at the next critical value a„+i. Very 
rough estimates give the location of restoration of real eigenvalues in the interval 
(n7r/2a, (2n + l)7r/4a). 
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(a) o = (b) a = 3.1 

Figure 4: Graphs of l{k) (full line) and r(fc) (dashed line), a = n/A, (3 = —0.5. 



In view of the presented arguments, only one complex conjugated pair can 
appear in the spectrum for fixed a and /3 in the neighbourhood of 0^+0^, and 
the other eigenvalues approach (n7r/2a)^ as the distance from + increases. 
Moreover, for fixed /3, the enlarging of a results into the shift of eigenvalues from 
almost Neumann ones (n7r/2a)^,n g N, to Dirichlet ones ((n + l)7r/2a)^, n € N, 
for a large. 

Finally, the equation for eigenvalues and eigenfunctions are found in a stan- 
dard way. The general solution Acos{kx) + Bsm{kx) of —-0" = fc^fA is inserted 



into boundary conditions (4.17) and the condition for existence of non-trivial 
solutions A,B is the eigenvalue equation (4.19). □ 



Figures [5j |6j represent the a-dependence of the first four eigenvalues as ob- 
tained by a numerical analysis of (4.19). The numerical results confirm the 



above described behaviour. Let us remark that if /3 is positive, then the graph 
of r{k) is reflected by the x-axis and the effect of loosing intersections is not 
possible, hence the spectrum remains real. 



4.4.2 Positive curvature 

The eigenvalue problem for the Hamiltonian reads 



ix) -\- ianx^' ix) ^ — tjj{x) = k^tplx) in (—a, a), 

cos-^ X 

V''(±a) + ia-0(±a) = 0. 



(4.26) 



Solutions of (4.26) can be written down in terms of associated Legendre func- 
tions Pi^',Q" 

V'(x) = CMx) -f C2i'{x) = CiPi'")(sina;) + C2Ql'"^(sinx), (4.27) 



where 



V := - (\/l + 4A-l 



(4.28) 
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Figure 5: a-dependence of eigenvalues, zero curvature, a — 7r/4, f3 — 0.5. Red, 
green, and blue colour corresponds to m = 0, 1, 2 respectively. 




C2{aM-a) - iV-aC-a)) = Ci(-aV'i(-a) + i^'ii-a)). (4.29) 



Inserting the general solution (4.27) into boundary conditions in (4.261 and 
consequent search for non-trivial constants Ci , C2 yields the eigenvalue equation 



V'i (a) + iatpi (a) ip^ (a) + \a1p2 [a] 
•4j'i{—a)+ia%lJi{—a) ijj'2(a) + iaijj2{ci) 



= 0. 



(4.30) 



In order to analyse the spectrum in more details, we transform the Hamil- 
tonian iJ^J^i) into a unitarily equivalent operator of a more convenient form. 
The proof of the lemma is a straightforward calculation. 

Lemma 4.3. The unitary mapping : L^{J2,'^x) — L^( J2, di^(+i)) 

(C/(+i)V)(a;):=(cosa:)-5V(a:) (4.31) 
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transforms H^^^^-^ (a, 0) to 

t^r+i)^iT+i)("' = ^i°(o)("' 5 tana) + 

^(+i)(^) 



8m^ — 3 — cos 2a; 



(4.32) 



(4.33) 



Equipped with the equivalent form of the Hamiltonian, we prove the follow- 
ing result. 

Proposition 4.4. For every m G Z there exists a real number ^^i^ such 



that all eigenvalues X with Re A > o.^^ real and simple (i.e. the 

multiplicity being one). The eigenvalues with Re A < ^"^i-^ can be complex, 
ordered in complex conjugated pairs. 



Eigenvalues are determined by equation (4.30) and eigenf unctions can be 



written in the form (4.271 with (4.291. 



Proof. Let us consider the transformed Hamiltonian (4.32) and forget about the 
potential for a moment, i.e., we understand the potential V^^^^ as a perturbation 
of H^gJa, ^tana). Since tana is positive under the assumption a < tt/2, 
the reality of the spectrum is guaranteed by Proposition |4.3| 1. The potential 
represents a bounded perturbation and it can shift eigenvalues only by C*]! V^^™]^^ || . 
Here the constant C comes from the estimate of the norm of the resolvent 



|i?i%)(A)|| 



< 



C 



which is valid for H^Q^{a, |tana) due to the similarity to a normal operator 



(c/. Theorem 4.1). The separation distance |A„+i — A„| of eigenvalues (ordered 
with respect to the real part) of the unperturbed operator _ff^Qj(a, | tana) 
grows to infinity and two eigenvalues must collide at first to create a complex 
conjugate pair. Hence, the perturbed operator cannot have more than finitely 



many complex eigenvalues. Recall that due to T'T-symmetry (Corollary 3.2 1 
the complex eigenvalues come in complex conjugated pairs. □ 

Remark 4.6. In other words, we detected the effects of positive curvature. It acts 
as the adding of real bounded potential V^^^ and real "/3 like" term in the 
boundary conditions to the zero curvature Hamiltonian -ffj^o)- The positive 
^ tana term is decisive for the behaviour of the spectrum, the bounded poten- 
tial V^+i) can affect substantially only the lowest eigenvalues. Nonetheless, we 
conjecture that the spectrum remain real for every m G Z. 

A numerical analysis of the equation ( |4.30 ) for A = is presented in Fig- 
ure [7] Obvious similarity with Figure [5] supports the perturbative results. 
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Figure 7: a-dependence of eigenvalues, positive curvature, a — it/4. Red, green, 
and blue colour corresponds to m = 0, 1, 2 respectively. 



4.4.3 Negative curvature 

The eigenvalue problem for the Hamiltonian H^_^-^ reads 



-4>" {x) — ianhxip'{x) + 



-^i;{x) = k i}{x) in (-a, a), 

cosh x (,4.C)4j 



V''(±a) + iaV'(±a) = 0. 



The solutions of (4.34) can be again expressed via associated Legendrc functions 
pI^\ Qu^\ but they have a little bit more complicated form then ( |4.27 ): 



where 



/i := im 



V cosh X 
1 



Qif^(tanha;) 
V cosh X 



Vl - 4A. 



(4.35) 



(4.36) 



Relations between Ci,C2 can be obtained from equation (4.29), however, with 
ipi,'4'2 corresponding to the negative curvature solutions (4.35); the same is true 
for the eigenvalue equation (4.30). 

To explain the behaviour of the spectrum in a deeper way, we use the same 
strategy as in the positive curvature case. The eigenvalue problem (4.34) can 
be transformed by an analogous unitary transformation leading to a modified 
zero curvature eigenvalue problem. 

Lemma 4.4. The unitary mapping : L^(J2,da;) — > J2, dz/(_i)) 

(J7(_i)?/') (x) (cosha;)"3 ip^x) (4.37) 
transforms H^_T^^(pi,0) to 



Vd^iT-D = ^1%) - o tanh a) + , (4.38) 
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(a) Real part of A 



(b) Imaginary part of A 



Figure 8: a-dependence of eigenvalues, negative curvature, a = 7r/4. Red, green, 
and blue colour corresponds to m = 0, 1, 2 respectively. See animation! for an 
animated visualization of the a-dependence of the eigenvalues. 



where 



^("!i)(^) := 



8to2 + 3 + cosh 2x 
8 cosh^ X 



(4.39) 



Proposition 4.5. For every m G Z there exists a real number A7^,^ such that 



all eigenvalues A with Re A > -^^(^1) either real and simple (i.e. the algebraic 
multiplicity being one), or there is one complex conjugated pair of eigenvalues 
with real part located in the neighbourhood of + . The eigenvalues with 
Re A < A^'^-^-j can he complex, ordered in complex conjugated pairs. 



Eigenvalues are determined by equation (4.30) with'ipi,^2 from (4.35). Eigen 



functions can be written in the form (4.35) with constants Ci,C2 satisfying 



(4.29) with'ipi,ip2 from (4.35) 



Proof. The proo f is the same as in the positive curvature case, cf. the proof of 
Proposition 4.4 The unperturbed Hamiltonian i?^Q^(Q!, — ^ tanha) corresponds 
to the case analysed in Proposition 4.3 2. □ 



Remark 4.7. The curvature effect is now represented by the bounded real po- 
tential and the extra negative term — ^ tanha in the boundary conditions. 

A result of the numerical analysis of the eigenvalue problem is presented in 
Figure [8j The resemblance to zero curvature case with negative /3 in boundary 
conditions is obvious. 



4.5 Connected boundary conditions 

The connected boundary conditions are, by their nature, more complicated than 
the separated ones and moreover, they are given by three real parameters b, c, (j). 
Like for the separated boundary conditions, we can use the unitary transforma- 



tions J/(±i) introduced in (4.31), (4.37) to transform the problems to the zero 
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curvature case, however, with modified boundary conditions and with addi- 
tional bounded real potentials defined in (4.33), (4.39). The modification 



of boundary conditions is presented in appropriate subsections below. 

The spectrum is not analytically described so far even for the zero curvature 
model and it is beyond the scope of this article to proceed with this analysis. 
The main aim of this section is to show the effect of curvature, i.e., the trans- 
formation of curved models to the zero curvature case. Furthermore, we present 
some results of a numerical analysis for the 'lowest' eigenvalues: ^-dependence 
for selected values of b,c. It is important to note that, unlike in the separated 
case, we do not start with our parameters b, c, (j) from a self-adjoint operator for 
= 0, as it was the case for a = in the case of separated boundary condi- 
tions. We remark that the case b — c = 0, (f> — ±7r/2 corresponds to irregular 
boundary conditions and the spectrum of such operators is completely different 
from the cases presented here (c/. [52j). 

4.5.1 Zero curvature 



We impose connected boundary conditions ( 2.13p ) on the solutions of eigenvalue 
problem for H^if^^ {b, c, </>) and we obtain the following equation for eigenvalues 



- 2fc + 2fc cos(2aA:)Vl + be cos (p + {bk^ - c) sin(2afc) = (4.40) 
and eigenfunctions 

ip{x) = Ci cos{kx) + C2 sin{kx), (4-41) 
where the constants are further restricted by 



C2 ((-1 + Vl + fcce'"^) cos(afc) + bksm{ak)^ 

= Ci ((1 + Vl + &ce'^) sin(afc) - bkcos{ak)J . (4.42) 
Proposition 4.6. Eigenvalues X ~ k^ of H^^.^^Ab^c^cj)) are determined by equa- 



tion (4.401, eigenfunctions read (4.41) with (4.42) 



Eigenvalues for m ^ can be obtained by the shift A H- A + while the 
corresponding eigenfunctions remain the same. 

Figure [9] illustrates the behaviour of eigenvalues for a certain choice of pa- 



rameters. 



4.5.2 Positive curvature 

The solutions of the eigenvalue problem for H^^j^^^ {b, c, 0) with connected bound- 
ary conditions ( 2.13p ) are the same as (4.27) except the constants Ci,C2 now 
satisfy 

C2 (yTTb^e''f'M-a) - Ma) + bM~a)) 

= Ci (yrT6ce''^Vi i-a) - A (a) + bM-a)) ■ (4.43) 
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(a) Real part of A 



(b) Imaginary part of A 



Figure 9: 0-dependence of eigenvalues, zero curvature, a = 7r/4, 6 = c = 0.01. 
Red, green, and blue colour corresponds to to = 0, 1, 2 respectively. 



0.15 0.20 



(a) Real part of A 




(b) Imaginary part of A 



Figure 10; (/)-dependence of eigenvalues, positive curvature, a = 7r/4, 
b = c — 0.01. Red, green, and blue colour corresponds to to = 0, 1, 2 respec- 
tively. 



The equation for eigenvalues reads 



-Vl + b?.e''t'M-a) + Ma) - bip[i-a) -^/T+T^ e'*V2(-a) + Ma) - bM-a) 
-ctpii-a) - VlTbce-''f'i}[{-^a) + ij[{a) -ci/^al-a) - \/rT6^ e-'*i/'2(-a) + Ma) 

(4.44) 



Figure 10 illustrates the behaviour of eigenvalues for a certain choice of the 
parameters. 

We employ the unitary transformation C/(+i) introduced in Lemma 4.3 to 
map H^^^^^(b, c, (j)) to a zero curvature Hamiltonian. 



Proposition 4.7. The unitary mapping defined in (4.311 transforms the 

Hamiltonian i/j^^-i) i^, c, 0) to 



= 0. 



(4.45) 
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(a) Real part of A 



(b) Imaginary part of A 



Figure 11: ^-dependence of eigenvalues, negative curvature, a = 7r/4, 
b = c = 0.01. Red, green, and blue colour corresponds to to — 0, 1, 2 respec- 
tively. 



where V^™]^) is defined in (4.331 anrf iJj'^jjQj 
of ijj e W'^'^{J2) satisfying 



dx' 



with the domain consisting 



*(a) 
^(+1) 



B(+i)*(-a), with ^'{x) := 



y^l + bc) e'"^ - ifetana 




and 



(4.46) 



c — \/l + be tan a cos < 



jfotan^a y/l + bct 



ifotana 



Eigenvalues X — k of H^^_^^-^{b,e,(j)) are determined by equation (4.44) 
eigenf unctions read (4.27) with constants Ci,C2 given by (4.43). 



Remark 4.8. The boundary conditions (4.46) are T^T-symmetric, but they are 



no more 7^-pseudo-Hermitian. This result shows that although we reduced the 
problem to the zero curvature case (in the sense of previous sections), the in- 
vestigation of spectrum must be done with more general boundary conditions 
than T'T-symmetric and 7^-pseudo-Hermitian at the same time. 



4.5.3 Negative curvature 

The solutions of the eigenvalue problem for Hl^_^s^ (6, c, (p) with connected bound- 
ary conditions ( 2.13p ) are the same as in the separated conditions case (4.35), 
but the relation between constants Ci,C2 is given by (4.43) with 4'iy4'2 corre- 



sponding to the negative curvature solutions (4.35); the same is also valid for 



the eigenvalue equation (4.44). 



Figure 11 illustrates the behaviour of eigenvalues for a certain choice of 
parameters. 

Proposition 4.8. The unitary mapping ?7(_i) defined in (4.37) transforms the 
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Hamiltonian H^f^_^^(b, c, (j)) to 

d2 



where V^^-^{x) is defined in (4.391 and ^^n(o) • 
sisting of %l: £ VF^'^(J2) satisfying 



(4.47) 

iwii/i i/ie domain con- 



*(a) 



-B(-i)^(-a), with 'I'(x) 




and 



(4.48) 



^\ + 5c) e'"^ + i6tanh( 



c + Vl + fectanhacos (/) + jfetanh^ a \'^Y+~bce "^ + i6tanha 



Eigenvalues X — k of H"^_^^{b, (j)) are determined by equation (4.44 1 with 
'ipi,'ip2 from (4.35). The eigenf unctions read (4.35), where constants Ci,C2 are 
given by (4.43) with ipi,^2 from (4.35). 



Remark 4.9. The boundary conditions (4.48) are T^T-symmetric, however not 
■P-pseudo-Hermitian, as for the positive curvature case. Thus again, it is nec- 
essary to investigate more general boundary conditions in zero curvature eigen- 
value problem. 



5 Concluding remarks 

The goal of this paper was to introduce a new class of T^T-symmetric Hamilto- 
nians defined in curved manifolds and describe the effects of curvature on the 
spectrum. Although we were able to find these effects for both separated and 
connected boundary conditions, the absence of results on reality of the spectrum 
for the latter (even in the case of zero curvature) did not allow us to present the 
conclusions in an entirely descriptive and explicit way. Let us therefore sum- 
marize the main features of the model for the separated Robin type boundary 



conditions (1.3) only. 

In Table [2| we schematically (and very roughly) describe qualitative prop- 
erties of the spectrum we observed in the constant-curvature cases. The entry 
describing the positive curvature case includes our conjecture (supported by 
numerical analysis) that all eigenvalues are real. 

One of the most instructive results in the paper are probably Lemmata |4.3| 
and |4.4[ which enable one to understand the effect of curvature in terms of 
an additional effective potential and boundary-coupling interaction. For the 



s-wave modes {i.e. m = in the decomposition (4.7)) and infinitesimally thin 
strips {i.e. a <^ I), it follows from the lemmata that the positive and negative 
curvature acts as an attractive and repulsive interaction, respectively. This is in 
agreement with a spectral analysis of similar models in the self-adjoint case of 
Dirichlet boundary conditions [SSHSj. However, the additional boundary inter- 



action is not negligible for positive widths a, and its effect is actually completely 
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curvature 


spectrum 


eigenvalues 


zero 


M 


only some a-dependent, crossings 


positive 


M 


all a-dependent, no crossings 


negative 


C 


all a-dependent, crossings, 

creation and annihilation of complex pairs 



Table 2: A heuristic summary of our analytical and numerical analysis. 



opposite (c/. Remarks 4.6 4.7): the positive and negative curvature gives rise to 



an attractive and repulsive Robin-type boundary condition, respectively. The 
interplay between these two effects is further complicated by the presence of the 
repulsive centrifugal term for |m| > 1, and the numerical analysis confirms that 
the overall picture of the spectrum can be quite complex. 

It follows from previous comments and remarks that there remain several 
open problems, e.g. the proof of the reality of all eigenvalues in the positive 
curvature model. Nonetheless, we would like to mention also some other inter- 
esting directions of potential future research: the spectral effect of curvature in 
non-constant curvature and non-constant boundary-coupling functions setting, 
the existence of Riesz basis for such setting or models defined on unbounded 
domains (waveguides) in curved spaces. The last case can be viewed as a natural 
continuation of [TO] where a planar "PT-symmetric waveguide was studied. 



Acknowledgement 

P.S. is thankful to B. Mityagin for very valuable discussions. The work was 
partially supported by the Czech Ministry of Education, Youth and Sports 
within the project LC06002. P.S. appreciates also the support of CTU grant 
No. CTU0910114. 



References 

[1] Adams, R. A. Sobolev spaces. Academic Press, New York, 1975. 

[2] Albeverio, S., Fei, S. M., and Kurasov, p. Point Interactions: PT-Hermiticity and 
Reality of the Spectrum. Letters in Mathematical Physics 59 (2002), 227-242. 

[3] Albeverio, S., Gunther, U., and Kuzhel, S. J-self-adjoint operators with C- 
symmetries: an extension theory approach. Journal of Physics A: Mathematical and 
Theoretical 42, 10 (2009), 105205 (22pp). 

[4] Albeverio, S., and Kuzhel, S. One-dimensional Schrodinger operators with V- 
symmetric zero-range potentials. Journal of Physics A: Mathematical and General 38, 
22 (2005), 4975-4988. 

[5] Andrianov, a. a.. Bender, C. M., Jones, H. F., Smilga, A., and Znojil, M., Eds. 
Quantum Physics with Non-Hermitian Operators (2009), vol. 5, SIGMA. 



34 



[6] Bendali, a., and Lemrabet, K. The effect of a thin coating on the scattering of a 
time-harmonic wave for the Helmholtz equation. SIAM Journal on Applied Mathematics 
56, 6 (1996), 1664-1693. 

[7] Bender, C. M. Making sense of non-Hermitian Hamiltonians. Reports on Progress in 
Physics 70, 6 (2007), 947-1018. 

[8] Bender, C. M., and Boettcher, S. Real Spectra in Non-Hcrmitian Hamiltonians 
Having PT Symmetry. Physical Review Letters 80 (1998), 5243-5246. 

[9] Bender, C. M., Brody, D. C, and Jones, H. F. Complex Extension of Quantum 

Mechanics. Physical Review Letters 89 (2002), 270401 (4pp). 

[10] BORISOV, D., AND Krejciri'k, D. PT-symmetric waveguides. Integral Equations Operator 
Theory 62, 4 (2008), 489-515. 

[11] BouCHlTTE, G., AND PETIT, R. On the concepts of a perfectly conducting material and 
of a perfectly conducting and infinitely thin screen. Radio Science 24 (1989), 13-26. 

[12] BOULTON, L., Levitin, M., AND Marletta, M. a PT-symmetric periodic problem with 
boundary and interior singularities. Journal of Differential Equations (2009). 

[13] Chugunova, M., and Pelinovsky, D. Spectrum of a non-self-adjoint operator associated 
with the periodic heat equation. Journal of Mathematical Analysis and Applications 342, 
2 (2008), 970-988. 

[14] Clark, I. J., and Bracken, A. J. Effective potentials of quantum strip waveguides and 
their dependence upon torsion. Journal of Physics A: Mathematical and General 29, 2 
(1996), 339-348. 

[15] Davies, E. B. Spectral theory and differential operators. Cambridge University Press, 
1995. 

[16] Davies, E. B. An indefinite convection-diffusion operator. LMS Journal of Computation 
and Mathematics 10 (2007), 288-306. 

[17] Davies, E. B. Linear operators and their spectra. Cambridge University Press, 2007. 

[18] Engquist, B., and NeDelec, J.-C. Effective boundary conditions for electromagnetic 

scattering in thin layers. Rapport interne CMAP 278, 1993. 

[19] Evans, L. C. Partial Differential Equations. AMS, Providence, 1998. 

[20] Fring, a., Jones, H., and Znojil, M., Eds. Pseudo-Hermitian Hamiltonians in Quan- 
tum Physics VI (2008), vol. 41, Journal of Physics A: Mathematical and Theoretical. 

[21] Garcia, S. R. The Eigenstructure of Complex Symmetric Operators. Proceedings of the 
Sixteenth International Conference on Operator Theory and Applications (IWOTA 16) 
179 (2008), 169-183. 

[22] GiLBARG, D., and Trudinger, N. S. Elliptic Partial Differential Equations of Second 

Order. Springer- Vcrlag, Berlin, 1983. 

[23] Gray, A. Tubes. Addison-Wesley Pubhshing Company, New York, 1990. 

[24] Hartman, p. Geodesic parallel coordinates in the large. American Journal of Mathe- 
matics 86 (1964), 705-727. 

[25] Hebey, E. Nonlinear Analysis on Manifolds: Sobolev Spaces and Inequalities. AMS, 
New York, 2000. 

[26] Jacob, B., Trunk, C, and Winklmeier, M. Analyticity and Riesz basis property of 
semigroups associated to damped vibrations. Journal of Evolution Equations 8 (2008), 
263-281. 

[27] Jain, S. R., and Ahmed, Z., Eds. Non-Hermitian Hamiltonians in Quantum Physics 
(2009), vol. 73, Pramana - Journal Of Physics. 

[28] Kaiser, H.-C, Neidhardt, H., and Rehberg, J. Density and current of a dissipative 
Schrodinger operator. Journal of Mathematical Physics 43, 11 (2002), 5325-5350. 



35 



[29] Kaiser, H.-C, Neidhardt, H., and Rehberg, J. Macroscopic current induced bound- 
ary conditions for Schrodinger-typc operators. Integral Equations and Operator Theory 
45, 1 (2003), 39-63. 

[30] Kaiser, H.-C, Neidhardt, H., and Rehberg, J. On 1-dimensional dissipative 
Schrodinger-type operators their dilations and eigenfunction expansions. Mathematis- 

che Nachrichten 252 (2003), 51-69. 

[31] Kato, T. Perturbation theory for linear operators. Springer- Verlag, 1966. 

[32] Klaiman, S., Gunther, U., and Moiseyev, N. Visualization of Branch Points in VT- 

Symmctric Waveguides. Physical Review Letters 101, 8 (2008), 080402 (4pp). 

[33] Krejciri'k, D. Quantum strips on surfaces. Journal of Geometry and Physics 45 (2003), 
203-217. 

[34] Krejcirik, D. Hardy inequalities in strips on ruled surfaces. Journal of Inequalities and 
Applications 2006 (2006), (lOpp). 

[35] Kre.tciri'k, D. Calculation of the metric in the Hilbert space of a PT-symmetric model 
via the spectral theorem. Journal of Physics A: Mathematical and Theoretical 4I, 24 
(2008), 244012 (6pp). 

[36] Krejciri'k, D., Bi'la, H., and Znojil, M. Closed formula for the metric in the Hilbert 
space of a PT-symmetric model. Journal of Physics A: Mathematical and General 39, 
32 (2006), 10143-10153. 

[37] Kre.tciri'k, D., and Tater, M. Non-Hermitian spectral effects in a PT-symmetric 
waveguide. Journal of Physics A: Mathematical and Theoretical 4I, 24 (2008), 244013 
(14pp). 

[38] Langer, H., and Tretter, C. A Krein Space Approach to PT-symmetry. Czechoslovak 
Journal of Physics 54 (2004), 1113-1120. 

[39] MiKHAjLOV, V. Riesz bases in £2(0, 1). Sov. Math., Dokl, translation from Dokl. Akad. 
Nauk SSSR II4, 981-984 (1962) 3 (1962), 851-855. 

[40] Mitchell, K. A. Gauge fields and extrapotentials in constrained quantum systems. 
Physical Review A 63, 4 (2001), 042112 (20pp). 

[41] Mostafazadeh, A. Pseudo-Hermiticity versus PT-symmetry I: The necessary condition 
for the reality of the spectrum of a non-Hermitian Hamiltonian. Journal of Mathematical 
Physics 43, 1 (2002), 205-214. 

[42] Mostafazadeh, A. Pseudo-Hermiticity versus PT-syinmetry II: A complete character- 
ization of non-Hermitian Hamiltonians with a real spectrum. Journal of Mathematical 
Physics 43, 5 (2002), 2814-2816. 

[43] Mostafazadeh, A. Pseudo-Hermiticity versus PT-symmetry HI: Equivalence of pseudo- 
Hermiticity and the presence of antilinear symmetries. Journal of Mathematical Physics 
43, 8 (2002), 3944-3951. 

[44] Mostafazadeh, A. Pseudo-Hermitian Quantum Mechanics. arXiv:0810.5643v2, 2008. 

[45] Naimark, M. Linear differential operators. Part I: Elementary theory of linear differ- 
ential operator. Translated by E.R. Dawson. English translation edited by W.N. Everitt. 
New York; Frederick Ungar Publishing Co. XIII , 1967. 

[46] Naimark, M. Linear differential operators. Part II: Linear differential operators in 
Hilbert space. Translated by E. R. Dawson. English translation edited by W. N. Everitt. 
New York: Frederick Ungar Publishing Co. XV , 1968. 

[47] Rubinstein, .1., Sternberg, P., and Ma, Q. Bifurcation Diagram and Pattern For- 
mation of Phase Slip Centers in Superconducting Wires Driven with Electric Currents. 
Physical Review Letters 99, 16 (2007), 167003. 

[48] Rubinstein, J., Sternberg, P., and Zumbrun, K. The Resistive State in a Supercon- 
ducting Wire: Bifurcation from the Normal State. Archive for Rational Mechanics and 
Analysis 195 (2010), 117-158. 



36 



[49] RusCHHAUPT, A., Delgado, F., and Muga, J. G. Physical realization of -symmetric 
potential scattering in a planar slab waveguide. Journal of Physics A: Mathematical and 
General 38, 9 (2005), L171-L176. 

[50] SCHOLTZ, F. G., Gbybr, H. B., and Hahne, F. J. W. Quasi-Hermitian operators in 
quantum mechanics and the variational principle. Annals of Physics 213 (1992), 74-101. 

[51] SlEGL, P. Quasi-Hermitian Models. Master's thesis. Faculty of Nuclear Sciences 
and Physical Engineering, CTU Prague, http://ssmf.fjfi.cvut.cz/2008/siegLthesis.pdf, 
2007/2008. 

[52] SiEGL, p. Surprising spectra of PT-symmetric point intera<;tions. aj:Xiv:0906.0226, June 
2009. 

[53] SlEGL, p. The non-equivalence of pseudo-Hermiticity and presence of antilinear symme- 
try. Pramana - Journal of Physics 73, 2 (2009), 279-287. 

[54] Spivak, M. a Comprehensive Introduction to Differential Geometry, vol. IV. Publish or 
Perish, Boston, Mass., 1975. 

[55] Wachsmuth, J., and Teufel, S. Effective Hamiltonians for Constrained Quantum 
Systems. arXiv:0907.0351, 2009. 

[56] Weir, J. An indefinite con vcct ion-diffusion operator with real spectrum. Applied Math- 
ematics Letters 22, 2 (2009), 280-283. 



37 



